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Abstract: With both conventional and modern techniques the 
author undertakes to determine the proper systematica of the lady 
ferns of the genus Athyrium found in North America north of Mexico. 
A taximetric comparison of ninety-nine characters by the computer 
program working on a non-statistical optimization measure of cluster 
analysis (Rubin. 1966) is used to test the validity of conclusions 
reached by other methods. A semi-objective graphic evaluation of 
their phylogenetic relationship is considered by the ground plan 
method (Wagner, 1961). Both analyses produce results essentially 
similar to those reached by biosystematic studies, namely that there 
are two distantly related lineages in Athyrium, one comprising the 
silvery athyrium and glade fern, the other the common and the 
alpine lady ferns. The glade fern is found to be a highly specialized 
member of its group, as is also the American alpine lady fern. Thus, 
there are four distinct species of Athyrium with the single species 
of lady fern consisting of three subspecies with geographically 
distinct ranges. 


I. INTRODUCTION 

The experimental investigation and description of patterns of variation in or¬ 
ganisms is termed taxonomy. Taxonomy is based on multiple correlations. Tradi¬ 
tionally, these multiple correlations have been perceived by the intuition method, a 
process fraught with the possibility of errors. The taxonomists mentally assess all 
descriptive information which becomes available for each specimen, in comparison 
with similar data from every other specimen, and finally they break up the many 
specimens into groups which form a hierarchy. This cognitive process of trial and 
error without any formal sets of rules, is for the most part, neither ordered nor 
rationally inferential and thus does not unequivocally indicate how the decisions 
have been or must be reached. 

In recent years, there have been numerous attempts, both to tackle the theore¬ 
tical foundation of taxonomy and to improve objectivity and reliability of methods 
employed. To the former case, it is a “numerical evaluation of the affinity or simi¬ 
larity between taxonomic units and the ordering of these units into taxa on the 
basis of their affinities” (Sneath and Sokal, 1962). In other words, it is the discri¬ 
mination of classes or, given a system of classes, the diagnosis or identification of 
individual taxon (Jardine and Sibson, 1971). We have our thinking process formu¬ 
lated into an explicit set of rules governing our methods of classification (Lagendre 
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& Vaillancourt, 1969). Vigorous logic of the exact sciences like mathematics is 
applied to this particular branch of biology (Rogers, Fleming & Estabrook, 1967). 
The authors intend to amplify and objectify traditional intuitive processes practiced 
by trained taxonomists and hope that these methods will produce stable and precise 
taxonomies based on objective criteria. To the latter case, methods of reducing 
masses of data into understandable and workable units, either manually (Anderson, 
1941, 1949) or with the aid of a high speed electronic computer (Sokal and Sneath, 
1963; Wirth, Eastabrook and Rogers, 1966) have been developed. 

The concept in numerical taxonomy is, in fact, simple and straightforward, but 
the methods involved are procedures of rather rigid and exact nature. It is a way 
of putting similar objects together into groups and dissimilar objects into different 
groups, with the number of groups and the descriptive or discriminant characteris¬ 
tics of the groups initially unknown. The number of characters measured on each 
object should be large, and ideally they should be uncorrelated. There are now 
many diverse methods of measuring similarity and of analysing matrices of simi¬ 
larity values (Blackith Reyment, 1971). In order to choose among these methods 
one needs to know their relative merits and the kinds of (systematic) information 
they produce. Unfortunately, there is no absolute test to decide which method is 
best, in whatever way a systematist chooses to define best. 

Although the techniques of numerical methods in taxonomy are still being 
developed, they have been used more and more often both in the construction of 
taxonomic systems (classification) and as an aid in biosystematic investigations 
(Lockhart and Liston, 1970). In this paper, an attempt has been made to apply 
computer-aided taximetric methods to detect and order multiple correlations in 
North American Athyrium, with the hope that it will help in the delimitation of 
species and infraspecific units. For purpose of comparison, a more or less objective 
and quantitative method has also been tried in the elucidation of evolutionary 
phylogeny of Athyrium (the method of Wagner Trees). 

I hope that these combined applications of phenetic and phylogenetic approaches 
to systematic problems will inspire gt eater confidence than those deduced from one¬ 
sided consideration of traditional methods. 

II. THE PROBLEM 

The monumental work of Butters (1917) on Athyrium in North America distin¬ 
guished three species of North American lady ferns. The western American lady 
fern appeared to him to be very similar to European material ( Athyrium filix-femina 
(L.) Roth). The eastern American plants were taken as two species, Athyrium 
angustum Presl., and A. asplenioides Desv. However, his studies did not resolve the 
taxonomic problems of Athyrium in North America north of Mexico to the satisfac¬ 
tion of all concerned. 

Some botanists agreed with him that there are two species of lady ferns in 
eastern North America and one species with two varieties in the West. At the 
other extreme, some authors proposed that all American materials should be con¬ 
sidered as a single species with many different forms (Shaver, 1954) or as geographic 
varieties of a single species (Weatherby, 1936). Most flora writers followed Butters 
(e.g. Small, 1938; Tryon et at., 1953; Wherry, 1961) while others upheld the idea 
that there is only one species (Fernald, 1946; Morton, 1963). 

There is continuing controversy about the validity of the characteristics dis- 
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tinguishing the several taxa of Athyrium in North America. Some botanists, such 
as Wherry (1948), feel that they can differentiate the eastern and western Athyrium 
even in a cruising automobile. Others, such as Fernald (1946), have maintained 
that there is no means at all to separate the two eastern taxa of Athyrium correctly, 
even by using the microscopic characters as enumerated by Butters. 

In view of the fact that several successful studies of Athyrium in the old world 
have been worked out by using characters of pubescence (Holttum, 1947), spore 
(Sledge, 1962) and cytological information (Bir, 1965; Mehra & Verma, 1954) and 
the availability of newer taxonomic methods such a study seems practical since 
Athyrium (esp. the lady ferns) are widespread and conspicuous members of the 
fern flora which are within easy reach. It is therefore appropriate to re-examine 
such a common group of plants which cause so much difficulty to so many laymen 
and professional botanists alike. 

The present investigation is carried out with the above ideas in mind. We 
proceed with both conventional and modern principles and methods to attack the 
problem as will be exemplified by the paragraphs that follow. 

IU. THE PLANTS'” 

Lady ferns are pretty and decorative ferns belonging to the genus Athyrium 
(L.) Roth. They have been hailed as “the Queen of ferns” (Lowe, 1876) and numerous 
poets have devoted line after line of praise to this group of attractive and beautiful 
ferns (Clute, 1938). 

Generally, they are found in moist and shaded situations of all sorts, from deep 
woodlands to open swamps, stony pastures and sandy bogs. 

The rhizome of the lady fern is erect or decumbent, bearing broad and thin 
scales. Anatomically, the rhizome is a dictyostele with three to six or more meri- 
steles. In the majority of the taxa, the root trace departs before the opening of 
the leaf gap. The structure of the black and wiry roots is diarch. 

Athyrium has an amphiphloic siphonostele or solenostele with phloem on both 
sides of the xylem. The leaf trace diverges from the stele as a single strand. 
After transversing the cortex of the rhizome it enters the petiole. The leaf traces 
are usually strap-shaped in the petiole and these unite upwards in the rachis into 
a single U- or V-shaped vascular bundle. (See Fig. 1.) 

Athyrium are showy and large ferns. The herbaceous fronds are rather large, 
elegant and delicate. The blades are elliptical, lanceolate to oblong in shape and 
are bipinnate to tripinnate or more finely divided. Veins are simple or forking, 
with free tips. The stomates of Athyrium are of the attached type. There are 
three types of epidermal cells, angular, wavy and elongated-wavy. 

The scales of Athyrium are not clathrate and are usually linear, lanceolate, or 
oblong to deltoid ovate. Structurally, darkened or blackish scales found on the 

(1 ) The following taxa are included in the present study: (i) The European lady fern 
[Athyrium filixfemina (L.) Roth); (ii) The western American lady fern [A. filix 
femina (1..) Roth subsp. cyclosorum (Rupr.) C. Chr.); (iii) The northeastern American 
lady fern [A. filix femina (L.) Roth subsp. angustum (Willd.) Liew); (iv) The south¬ 
eastern American lady fern ( A. filix femina (L.) Roth subsp. asplenioides (Michx.) 
Liew); (v) The American alpine lady fern (A. distentifolium Tausch subsp. ameri- 
canum (Butters) Liew); (vi) The silvery athyrium (A. thelypterioides (Michx.) 
Desv.); (vii) The glade fern (A pycnocarpon (Spreng.) Tidestrom). 
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Fig. 2. A combined photomicrographic and diagrammatic presentation of the stipe anatomy 
of North American Athyrium. The leaf trace leaves the rhizome as a single vascular 
bundle (a), and is then transformed via (b) and (c) into two separate amphiphloic 
bundles (d) along the entire length of the stipe. Approximately at or above the first 
pinnae (pair) the bundles fuse to form a single U. or V-shaped bundle. Only xylem 
cells are sketched in the outlined vascular bundle (s). (f) and (g) are photomicro¬ 
graphs of (d) and (e) respectively, (h) shows the detailed structure of cells within 
one of the two bundles at (d). Photomicrographs were taken from free hand sections 
(appr. 60X). 
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basal portion of the stipe are normally fibrous, the yellowish-brown ones covering 
the young crozier are thin and m:mbraneous in texture. Some long slender and 
often blackish (sometimes green or straw-colored) cobwebby hairs are quite common, 
and stalked glands may also be found. 

The sori are short and oblong-lunate, curved, or horse-shoe-shaped, with a mem¬ 
braneous indusium of corresponding shape. Some sori are more or less round, with 
or without indusia. The indusium may form glandular-tipped multicellular hairs or 
eglandular hairs with or without long cilia. 

The sporangium of the lady fern is short stalked. Some of the stalks may 
be proliferating. The yellowish or brownish to blackish spores of Athyrium are 
kidney-shaped and monolete with a single linear laesura. They are perisporous, the 
exine varies from smooth, reticulate, scabrate to echinate. 

The gametophytes of Athyrium are of the usual greenish cordate type. Except 
for the glade fern, which possesses unicellular club-shaped hairs, all of our taxa 
have naked gametophytes. 

Athyrium are all found to have the same basic chromosome number of X equals 40. 

Very little is known about the chemistry oi Athyrium. But, they have been 
demonstrated to contain hydrocyanic acid, tannins, leuco-anthocyanidins and xanthones. 

To most general collectors, all three subspecies of American lady ferns are 
easily confused with the common wood ferns of the spinulose group ( Dryopteris 
spinulosa Watt., D. campyloptera Clarkson, D. dilatata Gray, D. intermedia Gray, and 
their hybrid complex). 

In D. spinulosa and D. campyloptera {D. dilatata Gray), the proximal superior 
pinnules are short as compared to the proximal inferior pinnules which are the 
longest In D. intermedia they are also of different lengths. Typically, the second 
proximal lower pinnules are shorter than the upper ones. All these differences 
are quite apparent, especially in the basal pinnae. Such is not the case for all three 
lady ferns, for they have both proximal upper and lower pinnules or segments equal 
in length. Further, these wood ferns are, in general, lacier and are sharper pointed, 
with almost spinulose teeth. 

In fertile materials, the wood ferns possess reniform indusia, whereas the lady 
ferns contain asplenioid, athyroid and diplazioid sori. In sterile materials, a frond 
helps solve the problem. They can be distinguished easily by the pattern of vascular 
strand of the stipe. The wood fern has numerous vascular bundles, typically five, 
along the whole stretch of the stipe, whereas the lady ferns contain just two in the 
lower part of the stipe, and these unite distally in the stipe or rachis giving rise to 
the first (pair of) pinna. 

Another source of confusion is the hay-scented fern, Dennstaedtia punctilobula 
Moore. (See Fig. 1. B) It so closely resembles the lady ferns in general appearance 
that, once, a Ph. D. candidate in taxonomy of vascular plants, collected a whole 
bunch of them for me when we were searching for lady ferns. However, the hay- 
scented fern bea"s copious gland-tipped hairs on the stipe and all over the plant. 
The lady ferns, on the other hand, are merely scaly to the naked eye. Moreover, 
the hay-scented fern is more lacy in appearance and has a cup-shaped indusium 
with the sporangia inside. 

IV. MATERIALS AND METHODS 

We employed one of the three measures of similarity most commonly used, the 
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similarity coefficients (Sneath, 1957; and modified by Rogers and Tanimoto, 1960) 
for our present study. The partitioning of a collection of OTUs (Operation Taxo¬ 
nomic Unit, i. e. any taxon or object to be classified) into various clusters in the 
“best” way proceeds by the optimization of a measure, a non-statistical solution to 
the problem. 

In essence, the method we used here can be summarized in the following sen¬ 
tences (Rubin, 1966, 1967). First, we define a function for evaluating (measuring) 
the amount of structure present in a given partition of a set of data (measured 
characteristics of each object), and then we try to find that partition for which the 
function is optimized. Thus, the data are organized into homologous or coherent 
groups for that particular partition. 

We call a criterion function to evaluate a set of groups a “splitting function”. 
The kind of splitting function we use is a geometric measure (not statistical) called 
the average object stability, and relies on similarity coefficients. Normally, we can 
perform a direct search for the optimum partition. Exhaustive search is feasible 
only for small object set. Here, a “hill-climbing” algorithm is used to shift through 
the partition in a sensible way, searching for the optimum partition by continuously 
trying to find a partition which is better than the best one previously found accord¬ 
ing to the measurement of stability. The algorithm does not require a numerical 
measure but only a decision as to which of two partitions is more highly structured, 
or more valuable. 

The splitting function is defined in terms of the “breaking value” of similarity 
(S*) and average inside and outside similarity values. The optimum partitions for 
different values of S* define a hierarchy which we call a “tree”. If S* is small 
(approaching zero), the number of groups formed will be few and large in size; if 
S* is large (approaching 1), the groups will be numerous and small. Thus, a den¬ 
drogram can be formed by using several different values of S*. 

The methods of experimental design, data accumulation, and data processing 
may be briefly listed as follows: 

1. Selection of OTU: 

The first step is to choose the specimens to be classified. Here, normal and 
mature materials of Atlnyrium from all over continental North America north of 
Mexico were used in this study. For checking purpose, materials from Eurasia 
were also chosen for the present study. Both herbarium specimens and living plants 
were used. 

2. Selection of Character: 

The characters of each specimen were recorded. In order to achieve reasonably 
stable and objective classification, ninety nine characters for each of the 170 speci¬ 
mens were measured. (See appendix I and II.) Only characters which vary over 
specimens studied are included in the character. The characters used in the present 
study were derived from various portions of the life cycle of Athyrium. Some suc'i 
as the length of the blade, were easily and accurately measured. Others, such as 
the size of the spore, required time-consuming procedures. These data are numerical 
representations of all the characters that have been considered of significance for 
the classification of the taxon involved. Most of the characters selected are those 
used by previous investigators. Several are new and are either too difficult to be 
routinely studied or have been unaccessible to previous authors. 
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Acquiring the necessary data on the 170 plants took more than half a year of 
intensive effort. Further time was required for the trial run of the program in the 
Computer Center, Columbia University and also for key-punching all these data 
onto standard 80 columns IBM cards. 


3. Estimation of Similarity Between OTIT’s: 


The program uses a Q-technique, in its first stage to compute a similarity matrix 
in which all the objects are compared with each other. Mixed discrete, rank-ordered, 
and continuous data are used. Two different similarity coefficients are chosen: 
Fractional Match and Tanimoto Match, where 

N u mber of Matche s__ _M 

Total Number of Characters N 


Fractional Match Coefficient, Su = 


Tanimoto Match Coefficient, Su = 

For continuous data (also rank-ordered data) a simple and practical way of 
handling match and mismatch is by computing a fractional match, /, and thus, the 
“number of matches” need not be an integer. Each continuous character is specified 
with two values: 

a. Lower limit of difference —if two objects differ within this limit, treat it 

as a match (i. e. add 1 to M). 

b. Upper limit of difference — if two objects differ beyond this limit, they are 

considered not to match (i.e. add nothing to M). 

All continuous data will be mapped into the interval (2, 255) by the transformation: 

X> = 128 + JV,.„ (2=^) 


Where U = mean of character 
a = standard deviation 

Nj,, = number of units into which one standard deviation gets 
mapped 

X = original value of character 
X' = transformed value of character 

Under normal conditions, Ni,, assumes the standard value of 16, with mean 128 and 
standard deviation 16. 

For characters with non-applicable states, modified similarity coefficients are 
used to compute similarity: 


Modified Fractional Match Coefficient = n 
where n, = number of matches 

m = number of mismatches. 

Accordingly, the Tanimoto Match Coefficient becomes: 


Modified Tanimoto Match Coefficient, Su = 


2i.it. -t 


>r non-applicable states characters are treated according to the following rules: 
If both objects have non-applicable states for a character treat neither as a 
match nor mismatch (i. e. add nothing to either n, or tu). 

. If one object has the non-applicable state, and the other has an applicable 
stat: or value, treat this as a mismatch (i.e. add 1 to nj and nothing to 


n,). 

This definition is identical and applicable with the previously given ones: 
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a. For discrete characters, if they match, add 1 to n„ if they do not match, 

b. For a continuous variable, treat as part match, part mismatch (i.e. first 

compute the fraction of match, /, then add / to n„ 1 —/ to n d ). 

The similarity function which is determined for every pair of specimen varies 
between zero and one. Two exactly similar OTU’s have a value of one. Two com¬ 
pletely dissimilar OTU’s have a value of zero. 

After calculating the similarity coefficient for each pair of specimens, the com¬ 
puter prints out this data in an OTU by OTU similarity matrix. It also computes 
the average similarity and provides this information. 

4. Cluster Analysis For Organizing Objects Into Groups: 

The main clustering technique used in the present study was developed by J. 
Rubin for IBM 360 (Rubin, 1955). The detailed system of computer program deve¬ 
loped for these methods is explained in Rubin (1967). The program is written in 
Fortran IV with a portion of it in Assembler language. It uses 128K byte memory 
of an IBM 360/75. Practical biological taxonomic application of this method can 
best be verified in Prance et al. (1969). 

The program attempts to partition the objects into groups and keeps changing 
the membership of the groups to maximize the splitting function called the “average 
object stability” by the “hill climbing” algorithm, thus arriving at the best possible 
way of grouping the specimens. 

The object stability is the term which describes a function relating the average 
similarity of an object for a group which it is a member and its maximum average 
similarity for any outside group including the empty group: 

Object Stability, O, = 

where M = average similarity of each object to each group 
(average inside similarity) 

M = maximum average similarity to any other group 
(maximum average outside similarity) 

S* = breaking coefficient of similarity. 

The first term at the right side of the equation is the normalized amount by which 
the inside attraction (average similarity between the object and the members of the 
group) exceeds the breaking value, and the second term is the normalized amount 
by which the breaking value exceeds the maximal outside attraction. The value is 
either zero or negative. This definition of object stability has been proven to be 
better than the simple relationship, Oi = M — M. 

Based on the object stability measured, each object can be labelled as stable or 
unstable — stable with a positive object stability and unstable with a negative object 
stability. This information is printed out in the output. 

The splitting function, which is defined as the average object stability, is com¬ 
puted over all objects in the data set. 

Average Object Stability = £ 0/ 

where n = total number of objects 
Oi = object stability 


Oct., 1972 Liew—Studies on North American Lady Ferns and their Allies 199 

We will split a set of objects into groups for a certain partition value if it satisfies 
the following two factors: 

a. the average similarity of the objects inside the group exceeds a certain 

breaking coefficient (S*), and 

b. the average similarity value between the members of different groups is as 

possible (less than, or at most, equal to S*). 

The splitting function is used to evaluate the degree of organization achieved by 
the grouping. It measures the value of this grouping and then decides which of 
every two groupings is better. The “hill climbing” procedures described below 
examines every possible way of grouping the objects, and selects the best one that 
optimizes the splitting function. 

With a given initial partition the effect of moving each object into every other 
group except the one it is in is calculated. If the new partition so created is better 
(i. e. gives a lower average similarity value between groups), the object is moved, 
otherwise, we move it back to the group it was in originally, i. e. we retain the 
better of the two partitions. We also try to move the object into the empty set, i.e. 
into a group by itself. Thus, a single object gets moved into the best possible 
group while other objects remain where they are. We then move on to the next 
object, and repeat the processes. If asked for, the effects of moving an object from 
one group to every other group (or only to the best group) will be printed out to 
enable one to figure out a partiton better than the one the program would reach. 

After processing every object, we start from the beginning again. When no 
movement of a single object can improve our "best” partition, we reach a “local 
maximum” which might be different from the true optimum partition. There may 
be more than one local maximum (the maximum allowable number of groupings is 
10 ). 

We use four types of heuristic techniques to overcome local maxima which are 
not optimum partitions: 

1. Forcing pass — By forcing each object out of the group that it is in, and 

into the outside group with highest average similarity attraction for it. 
The value of the partition created after all objects have been reassigned 
is then calculated and compared to the original partition. Forcing passes 
are repeated until they produce no improvement. 

2. Reassignment pass — Starting with the best known partition P, and reassigned 

each object to the group with the highest partition for it. The process is 
repeated until no improvement is noted. 

3. Forcing with suppression pass — These are identical to forcing passes, except 

that we do not permit moves into a group which was non-empty at the 
beginning of the forcing pass. (In effect, we are reclustering each group 
without taking into account the attraction from other groups.) 

4. Forcing without restoration pass — After forcing each group we do not re¬ 

store the best known partition but simply hill climb to a new local maxi¬ 
mum from whatever point is reached. If the local maximum is not an 
improvement over the previous ones, the procedure is terminated and the 
optimal partition is printed out. 

Each object is finally placed into a group where the final result is the greatest 
object stability that can be achieved by the partitioning. The final optimal partition 
is printed out together with the groups formed and the membership in each group. 
This is the result we are looking for. Here, the most stable object within a group, 
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i. e. the one with the highest object stability, can be considered as the most “typical’’ 
objective of the group. 

5. Computer Output: 

The computer is programmed to print out some important information in the 
following sequences: 

1. Similarity matrix—Fractional Match or Tanimoto Match. 

2. The number and membership of maximum — One or more local maxima 

reached and the groups so formed. 

3. Final set memberships for the best partioning value — This is taken as the 

only significant result from which taxonomic inferences are drawn. 

4. Effect of individual moves — Only the effect of the move into the “next best" 

set or all effects by moving each object into every group of the final 

optimal partition other than the one it is in. 

5. Average similarity between groups — The diagonal elements represent the 

“internal cohesion” of each group. It informs us which groups are well 

structured and which ill-formed. 

6. Average similarity of each object to each group —An object has a high 

average similarity to the group it is in, and low to all other groups. 

7. Description of stable elements — Object stability is calculated according to 

the object stability equation and each parameter involved is printed out. 

Each object in a group is labeled as stable or unstable. 

8. Variable analysis for each group —If asked for, a character distribution 
analysis for all three types of variables will be printed out. It includes frequency 
counts within each group for continuous data. 

All computations were carried out at the Computer Center, Columbia University, 
New York, with IBM 360 under release OS 18. The 170 objects data set was run at 
breaking coefficients (S*) from .30 to .80 with initial interval of .10. Finer divisions 
were used when it was considered necessary for producing meaningful groupings. 
Dump cluster or residue group is not allowed. Actual computing time was about 
three hours computer time. 

V. RESULTS 

We have two sets of results, one from the Tanimoto Match (T.M.) and the other 
from the Fractional Match (F.M.) of similarity coefficients. 

A. Tanimoto Match Coefficient (T.M.) 

Five levels of S*, .3, .4, .5, .6 and .7 are used to construct the “tree” for this 
measure. The resulting "tree” is depicted in Fig. 4. 

At .3, two groups are formed from the 170 specimens. One group includes A. 
pycnocarpon and A. thelypterioides; the other, the lady fern complex, includes both 
A. distentifolium and A. filix-femina sensu lato. The next level of S* produce 4 
groups, essentially by dichotomous branching of the two previously formed groups. 

At this level, A. pycnocarpon is separated from A. thelypterioides and A. disten- 
folium is also separated from A. filix-femina sensu lato. 

At S* = .50, seven groups are formed. A. pycnocarpon and A. thelypterioides 
continue to be separated; these two, in fact, retain their identity from S* = .40 on¬ 
wards to the highest S* value we used. The alpine lady fern differentiates into 
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two groups, one, the leafy alpine lady fern of Eurasia, the other, the leafy and 
skeleton-like specimens of North America. The lady fern complex splits into three 
groups, one containing materials from the northeastern American region, another 
from the southeastern American region, and a third large group containing materials 
from the western American region and also Eurasia. 



The Tanimoto Match coefficient is used in computing 
similarity of all OTUs is .413. The hexagons represcr 
tical line indicate connection for clusters formed al 
serted diagram is a simplified "t 


s of breaking coefficients (S*). 


differ! 


n Norl 


Athyrium pycnocarpon-, TH: A. thelypterioides ; 
Dl: A. distentifolium of North America; NE; 
A. filix-femina subsp. asplenioides-. SW: A. filix- 
tbsp. cyclosorum of North American southwest; NW: A. filix-femina subsp. 
» of North American northwest; EA: A. filix-femina of Eurasia. 


oximate. 

DI-EA; A. distentifolium of Eurasia 
A. filix-femina subsp. angustum-. SE 


At S* = .60, the program produced 9 groups. Six of the seven groups of the 
previously partited groups remained unchanged except the large group. It is 
separated into three smaller groups one containing specimens of Eurasia, one of 
northwestern America and the other southwestern America. 

This result seems provocative. Since they have enough morphological similarity, 
both northwestern and southwestern American lady fems are considered to be a 
single taxon despite their only difference in the spore character. This taxon has 
long indusial cilia which are eglandular. It has two kinds of spores (the north¬ 
western plant has smooth to papillate and yellowish spores only; the southwestern 
plants are mixed). When these characters are combined with the fact that the 
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western plant has gametophytes with close or approximate inner wing (as are other 
members of the American lady fern complex), the western American plant is distinct 
enough from the Eurasian materials to warrant it a subspecies status. 

The western American plants are more different from all eastern plants than 
from European materials. In this respect, the total absence of stalked glands, the 
long indusial cilia (eglandular), the short erect or suberect nature of the rhizomes 
etc. make them distinct from eastern American taxa. However, these characters, 
when taken singly, are generally not sufficient enough for subspecies differentiation, 
but when combined together, provide enough evidences for the proper identification 
most of the time. 

The results at S* = .50 and S* = .60 seem significant to the present study. Since 
our interest is mainly in the Athyrium of North America a simplified diagram using 
only North American materials is used to depict the dendrogram unequivocally. 
The average similarity between each group and also the intent homogenity of each 
group (diagonal values of the matrix) are given in Table I. For S* values over .60, 
the number of groups formed is too numerous (25), and memberships within each 
group are incomprehensible. They are therefore omitted from the diagram. 

The computer output for each final partition and memberships is given in the 
Appendix. Typical computer output containing useful information is given at S*=.60 
for both T. M. (9 groups) and F. M. (7 groups). 

Table I. The average similarity between each North American group of Athyrium 
at two values of breaking coefficients, S*, using Tanimoto Match for 
computing similarity coefficient. The diagonal values 
represents internal homogenity of each group. 
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B. Fractional Match Coefficient (F. M.) 

Five values are used as input of breaking coefficients, from .30 to .70 at an 
interval of .10. The output, drawn as a dendrogram, is given in Fig. 4. 


S 



Fig. 5. Formation of clusters for different values of breaking coefficient (S*). Fractional 
Match coefficient is used in computing similarity matrix. The average similarity of 
all OTUs is .524. Detailed information of the computer output on optimal partition 
with breaking coefficient at .60 (seven clusters) is given in the Appendix. Legends 
are as in Fig. 4. 

It is evident from the “tree ’ that at S*=.30 all of the 170 specimens were 
separated into two large groups, one comparising A. pycnocarpon and A. thelypterioides, 
and the other, A. distendfolium and the lady fern complex. At S*=.40 the grouping 
pattern does not change. This result manifested the fact that both groups are more 
similar within themselves than between each other. We will appreciate such a 
grouping better if we take into consideration that Butters (1917) grouped both A. 
pycnocarpon and A. thelypterioides under the same genus Diplazium. Anyone who 
has some acquaintance with lady ferns knows that morphologically common lady 
ferns and Alpine lady ferns look much alike in several aspects. 

The internal differentiation of the group of lady ferns starts at S*=.50, where 
A. distend folium splits out of the grouping. A. distendfolium remains a separate 
group all the way through higher levels of S*. The rest of the common lady 
ferns form a group by itself and continue through S*=.60. By that value of S*, A. 
thelypterioides and A. pycnocarpon separate from each other and each of them remains 
isolated throughout the rest of the S*. 
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At S*=.70, we have seven groups. Each of the group of A. pycnocarpon, A. 
thelypterioides and A. distentifolium is a separate group by itself. The largest group, 
comprising the common lady ferns, separates into four groups, which incidentally, 
correspond to the geographical distribution of the groupings. Thus, we have the 
northeastern, northwestern and southeastern lady ferns. 

The result is very interesting. In checking through the average similarity of 
each group in Table II, it is clear that northeastern, southeastern should be maintained 
as be separate groups. They have between group similarity values at or below .60. 
For northwestern and southwestern lady ferns, it is higher than .70. We thus 
conclude that the western materials should be grouped together as one group instead 
of two. This provides part of the basis for the decision on subspecies grouping 
within common lady ferns. 


Table II. Average similarity of each group at two values of breaking coefficients, 
S* (.60, .70). Fractional Match Coefficient is used in computing similarity. 
Notice that the diagonal values represent the internal 
cohesion of each group concerned. 

A.) At S« - .60. four groups. _ 

Taxon | M *^ r j StahMty j pycn0 [ ,hely | dist | FernB 

.71 .883 

.84 .560 .936 j 

.61 .281 281 .843 

.15 .295 .343 .424 .661 


pycno 15 2 

thely 15 j 4 

dist 20 1 

Lady Ferns 120 3 


B.) At S* = .70, seven groups. 


Taxon 

Membe 

Set 

Stability 

pycno 

h he ' y _ 

dist | ff-NE | ff-SE 

ff-NW 

ff-SW 

pycno 

15 

2 

.61 

.883 





thely 

15 

4 

.79 

.560 

.936 




dist 

; 20 

1 

.48 

.281 

.281 

243 



ff-NE 

30 

7 

.46 

.305 

225 

.356 .838 



ff-SE 

30 

3 

.28 

.334 

285 

.393 .587 .783 



ff-NW 

30 

6 

.51 j 

.253 

297 

.456 .595 .547 

.860 


ff-SW 

30 

5 

.58 

.290 

.365 

.491 | .577 .601 

1 - 713 

.879 


VI. DISCUSSION 


We performed two different kinds of similarity matrix calculations, utilizing the 
same procedure of clustering technique in both. The final optimal partition at each 
level is drawn to a hierachy "tree”. The interrelationships of each group are not 
shown diagramatically but are shown in Table I and Table II. 

All of the 170 OTUs are grouped nicely into clusters which agree reasonably 
well with groups drawn by close observation both in the field and in the herbarium. 
All members of each group are found to be stable within their respective groups, 
as measured by object stability. 
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The similarity value between each pair of OTU obtained in the Fractional Match 
is higher than in the Tanimoto Match. So is the average similarity coefficient, 
which in the former has a value of .524 and .413 in the latter. Accordingly, the S* 
value required to generate a certain number of groups should be lower in Tanimoto 
Match than in Fractional Match. This is exactly the case we encountered. Take 
for instance, at S* = .60, Tanimoto Match generate nine groups while Fractional Match 
only seven. Seven groups are produced in the Tanimoto Match at S* equals .50. 

Another minor difference is that each of the A. pycnocarpon-A. thelypterioides and 
A. distentifolium-A. filix femina groups branches dichotomously into its constituent 
members in the Tanimoto Match at S* = .40 while in the Fractional Match they 
diverged into two groups at S* = .50 and seven groups at .70. 

Both measures differentiate the Eurasian materials of A. distentifolium and A. 
filix-femina proper from the North American materials. Both measures also maintain 
that there are three basic groups in the Athyrium of North America. One group 
contains A. pycnocarpon and A. thelypterioides. The next one is A. distentifolium all 
by itself. A third group involves the whole of A. filix femina sensu lato. 

Both measures depict that A. distentifolium is closely related to A. filix femina 
sensu lato and rather distantly associated with A. pycnocarpon and A. thelypterioides. 
But A. distentifolium is separated from A. filix femina sensu lato at the first step of 
their further divergence. This conclusion has been reached independently by both 
measures. 

The main problem remains for the grouping of North American lady ferns 
complex. In T. M. measure the whole group splits into three groups at S* = .50. 
While both northeastern and southeastern taxa retain their identity at .60, the 
western and Eurasian lady ferns diverge into three groups, one Eurasian, one north¬ 
western, and another southwestern. Ignoring Eurasian materials, we have in North 
America four groups of lady ferns. Exactly the same number of groups is obtained 
with F. M. at S* = .70. 

It is therefore tempting to suggest that North American lady fern should be 
arranged into four natural groups. However, as stated in the result on F. M., the 
intergrouping similarity between northwestern and southwestern taxa is very high. 
Both are grouped together at S* = .50, and separated only at S* = .60 in T. M. 
Morphologically, it is not possible to separate the two except by using characteristics 
of spore sculpture. It is here proposed that both should be considered members of 
the same ultimate taxon. 

As can be seen from the dendrogram of F. M., both A. distentifolium and the 
common lady ferns occupy a long vertical stretch from at or below S* = .50 to at 
or above S* = .60 it is accordingly sensible to regard them as natural groups. Rubin 
(1967) wrote: “If each optimal partition over a long stretch i. e., over a relatively 
long interval of S* should contain a group which does not change from level to 
level then this group should be considered as a candidate for a natural group.” We 
may regard them as species. 

This is quite consistent with the investigation of Prance (1967) and Prance et 
al. (1969) in which he was using .40 and .50 as the branching coefficient for delimiting 
genera. Since our plants separate at S* values greater than .60 it is logical to 
recognize them as subspecies. The three subspecies thus delimited in North America 
are northeastern, southeastern and western lady forns. 

The introduction of numerical taxonomic technique here does not entail a com- 
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plete break through with the past. It does, however, represent a more refined and 
precise method in the detection of phenotypic clustering and discontinuity and hence 
the limits of taxa. 

The phylogenetic approach which follows this section generated results which 
correlate very well with the conclusions here drawn. Thus, this method of computer 
taxonomy does reflect, to a certain extent, phylogenetic relationships among Athytium. 
In this connection. Prance (1967) commented that “The technique is obviously a 
powerful method of clustering because of the way it has denoted most of the genera 
easily” and further in 1969, he wrote that “this method clearly confirmed the main 
thesis for making changes at generic level.” 

VII. PHYLOGENETIC CONSIDERATION 

Recently, numerous methods of objective estimation of phylogeny from similarity 
relationships among existing taxa have been developed by several investigators, 
both from morphological (Camin and Sokal, 1965; Wilson, 1965; Cacalli-Sforza and 
Edwards, 1967) and also from molecular biological standpoints (Fitch and Margoliash, 
1967, 1970; Horne, 1967). But, without enough paleontological evidence it is difficult 
to work out a phylogeny which truly reflects historical sequences of the past. We 
can only guess what has occurred in the past, based on information obtainable from 
existing materials. 

As has been mentioned before, most conventional taxonomic work on phylogenetic 
diagrams is based on taxonomic intuition and is, in general, quite subjective since 
only a few key characters, or a few characters that make for easy identification, are 
used for their constructions. Because key characters themselves are not always 
stable, more objective and repeatable methods based on large ensembles of data 
from all aspects of the taxon should be used. The numerical taxonomic method, 
employing a great abundance of characters, is essentially a development of the late 
fifties. It can be considered as a proper step toward reducing the subjectivity of 
systematic work to a minimum. The preceeding section is an attempt at this 
relative method in forming a classificatory hierarchy. 

For comparative purposes, a semi-objective graphic method for expressing evolu¬ 
tionary or phylogenetic deductions of relationship, the visual ground plan correlation 
method of Wagner (1961, 1969) based on the ground plan idea of Danser (1950) is 
used here. It is simple and based upon generally sound principles. This method is 
semi-objective since it is still the responsibility of the investigator to decide sub¬ 
jectively the relative degree of advancement or primitiveness of each character state, 
as well as the arbitrary values to assign for degrees of specialization. But uniform 
standards are used for all taxa concerned, and the construction of phylogeny is, 
therefore, not highly subjected to personal bias. It also utilizes a substantial number 
of characters and permits quantification of data. 

Numerous modifications of the original method have since been made (Whiffin 
and Bierner, 1972) and a computer program has recently been developed by Farris 
(1970) and have been applied successfully by Kluge and Farris (1969) on the evolution 
of anurans and the phylogeny of Gutierrezia (Solbrig, 1970). 

This method assumes three working hypotheses (Wagner, 1961): (1). Plants 

with a majority of similar characteristics have the same common ancestry. (2). 
Evolution proceeds normally in various directions, and different lines therefore 
change in different characters and different character complexes. (3). Evolution 
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occurs at different rates at various times and in different lines. 

The construction of a phylogenetic diagram based on these assumptions involves 
three broad phases: (1). The choice of characters — systematic or comparative 
analysis of the taxa in question to find out and understand their contrasting characters. 
(2). The determination of ground plans —to find out the character states common to 
all or most of the taxa in order to deduce the most probable ancestral or primitive 
states. (3). The diagramatic synthesis of phylogeny —a phylogenetic synthesis of 
data to assemble the taxa into a plausible graph according to their respective 
deviations from the basic ground plan and from each other. 

The working steps may be summarized as follows: 

1. To compare and study all the variable characters among all the taxa used; 

2. To determine the primitive condition by assuming that a character state 

which is present in the greatest number of taxa within a taxonomic group 
(and related taxonomic groups) represents the ground plan (generalized 
condition); 

3. To assign for each character the value of zero for the primitive condition, 

and one for the specialized condition, with intermediate states being given 
the value of 0.5; 

4. To list in tabular form the taxa and for each give the divergent values 

(specialization) from the ground plan, both for individual characters and 
in total; 

5. To determine the mutual character grouping between taxa and then arrange 

them in sequence according to these groupings on a concentric graph. 

The mutual character complexes (grouping) determine the radii and branchings, 
while the specialization (divergence) governs the distances. Taxa are connected to 
each other by their ensembles of common features, which are then plotted as the 
points of separation. Thus, the method displays a visual character grouping of the 
most probable common ancestors and outlines the pathways of phylogeny. 

In this study, 17 different characters were recorded. These are characters which 
differ in the different taxa and are believed to have evolutionary significance. 

An explanation (and scoring) of the primitive and advanced condition of these 
17 characters is given in Table III. The specialization indices for the taxa of 
Alhyrium in North America are tabulated in Table IV. Mutual character groupings 
are shown in Table V. 

The specialization index shows the relative degree of advancement for each 
taxon in one dimension. When this index is implemented with mutual character 
groupings (determined on the basis of the correlation of similar characters), a two- 
dimensional diagram portraying probable phylogeny of Athyrium can be constructed. 
This is shown in Fig. 5. 

In the target diagram two main evolutionary lines can be discerned. One 
involves A. pycnocarpon and A. thelypterioides, the other the A. distentifolium and A. 
filix-femina s. I. The latter consists of two major entities. At lower specialization 
index values (from 3.5 to 4.5) we have a group of A. filix-femina s. 1, and a very 
specialized shoot branching off from the ancestral line of Western American taxa 
of Athyrium. 

The degree of specialization in this group is rather interesting. There are only 
two major levels of specialization. One at specialization index equals 3.5 to 4.5 (5.5). 
Here are found the taxa of all the lady ferns, and at a slightly higher specialization 
index (5.5) the A. thelypterioides which arise from different ancestry. The other 
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Table III. List of characters and states used in the construction of phylogenetic 
relationship. Zero represents primitive condition, 

1 advanced and .5 intermediate. 


Habit of Rhizome and i 
Long creeping with 
Short erect with yt 
Shape of Frond Cluster 

Bunched compactly or in i 
Loosely bunched 
Complexity of Frond 
Leafy 

Skeleton-like 
Cutting of Lamina 

Bipinnatifid to quadripinm 
Once pinnate 
Cutting of Pinnae 

Simple 

Margin of Pinnae 

Notched and toothed 
Entire 

Forking of Veins 


Young Shoots 

oung shoots in front of the rhizome 
ig shoots in the center of rhizome 


K 


Vesture on stipe bases 
Mostly scaly 
Hairy and scaly 

Stalked glands (on rachis and costae) 
Absent 
Present 

Shape of Epidermal Cells 
Angular 

Elongate and wavy 
Elongate 


Shape of Sori 

Mostly aspienioid 

Mostly athylioid (and diplazioid) 

Roundish 
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Table III. Continued 


iracter and the description of the 
different states considered 


None 

Paraphyses (sori) 


Q 


Yellowish to brownish 
Sculpture of Perine 
Smooth to papillate 
Folded and with ridges 
Reticulate 

Perine wing of Spore 
Barely observable 


anastomosing 


level of specialization is 8. There are two groups at this level. One is A. pycno- 
carpott, the other, A. distentifolium. 

If the true phylogenetic relationships among North American species of Athyrium 
are correctly shown in Fig. 5 the following facts becomes evident: 

1. Athyrium of North America is a heterogenous group with two main evolution¬ 

ary lines. In one line, A. pycnocarpon is very specialized. The A. thelyp- 
terioides, though less specialized is still a little bit more specialized than 
all the lady ferns. It shares a common ancestry with A. pycnocarpon. 

2. A. distentifolium derived from the main evolving line of A. filix-femina s. 1., 

being more closely related to the western American taxa of Athyrium. 

3. This diagram separates the eastern American members of Athyrium filix- 

femina into two taxa quite distinctly, while the western taxa are closer 
in relationship to each other. They remain spearated only by the color 
of the spore (0) and their perine sculpture (P). 

These conclusions are very helpful in the study of the phylogeny of Athyrium 
in North America. Although such a diagram helps us to visualize a probable phy¬ 
logeny, it is not necessarily correct in all respects. Not all the pertinent characters 
could be used in the construction of the diagram. Many of them have their advanced 
condition not well defined. Many of the key characters could not be used simply 
because of the fact that they are rather subtle quantitative characters. Also, such 
a scheme does not easily accommodate reticulate and parallel evolution patterns if 
such patterns do exist. 

Moreover, the basic assumption of working steps number two (primitiveness) in 
the preceeding paragraph of this section is not always correct. If it were, we would 
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Table IV. Specialization indices for species and subspecies of 
Athyrium in North America**. 



species or subspecies. Mutual character complexes were then determined and utilized to 
plot on a concentric target diagram. 


Table V. Mutual character grouping used in the determination of radii and 
branching sequence of the phylogeny of North American Athyrium . 




Thely 

NW 

SW 

NE 

SE 

Dist 


al Character Grouping Specializ 

BDEFH1PQ 
BGIOPQ 
AJKLMO 
AJKLMP 
JKLMNO 
KLMNP 
ACKLMNPQ 


8.0 

5.5 


4.5 

3.5 


have to assume that apocarpy is more advanced among the angiosperms than syn- 
carpy or that gradate sori in ferns are more advanced than mixed sori, and 
eusporangiate ferns more advanced than leptosporangiate ones. 
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tion indices (distance) and mutual character groupings (radii and branching). Num¬ 
bered concentric lines represent levels of specialization. Capital letters refer to 
character states shown in Table VI and are given the first time the specialized condi¬ 
tion appears in an evolutionary line. Solid circles represent taxa, and empty circles 
represent hypothetical common ancestors. Note that the exact position of the connec¬ 
tion of the between taxa, the straightness of the lines, and the angles of divergences, 
are strictly diagrammatic and are not intended to denote a constant rate of divergent 
evolution among taxa. 

Despite all these drawbacks this method is an effort toward improving the 
objectivity of our determination of phylogenetic trends and relationships to the 
highest degree of probability. In the literature, this method and its modified versions 
have been successfully applied in several recent systematic investigations (e. g. 
Hardin, 1957; Hauke, 1963; litis, 1959; Keener, 1967; Kluge and Farris, 1969; Mickel, 
1962; Sokal, 1970; Stern, 1961.) 

As depicted in Fig. 5 the scheme produced by this method probably best repre¬ 
sents the most probable evolutionary relationships among North American Athyrium. 
Although the diagram cannot be claimed to represent the most accurate phylogeny 
of the taxa concerned, it will, however, serve as a useful first approximation, possibly 
it approximates the most parsimonious tree in the sense of Camin and Sokal (1965). 

It is worth pointing out that the phylogenetic relationships of Athyrium derived 
from this study resemble closely those derived from numerical taxonomic studies. 
It therefore validates directly the usefulness of both methods. 

The present author is in the process of constructing a phylogenetic tree of 
Athyrium in North America north of Mexico with the h lp of the computer program 
developed recently by Farris (1970). Due to difficulties in adopting the Forturan IV 
program for a CDC 3150 computer the computational procedure is not yet operational. 
The results will be reported later. 

VIII. COMPARATIVE DISCUSSION ON PHENETIC AND 
PHYLOGENETIC STUDIES 

Most authors believe that until all species of a genus have been studied any 
phylogenetic scheme proposed is speculative at best. Ideally, one should use all the 
data, from past to present, in his phylogenetic determination of evolutionary trends. 
However, in the virtual absence of any paleobotanical information in Athyrium this 
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cannot be done, and we have to utilize circumstantial evidences to postulate a ten¬ 
tative phylogeny. 

The phenetic and phylogenentic methods presented here are an effort to simulate 
the actual steps in the process of data manipulation for reaching taxonomic judge¬ 
ments. They help to make systematic conclusions more understandable, more reli¬ 
able and more objective. 

The most obvious characteristic of phylogenetic consideration which results in 
the production of target diagram is the emphatic way in which the relationships 
between groups are shown. This is not shown in the dendrograms of cluster ana¬ 
lysis, but is available in numerals in the tables of average similarity between groups. 

Both phenetic and phylogenetic studies indicate that there are two main groups 
of Athyrium in North America. They are the group of A. thelypterioides-A. pycno- 
carpon and the group of A. filix-femina complex with A. distnetifolium included in it. 

The group of A. thelypterioides-A. pycnocarpon forms a lineage which splits into 
individual members at a rather early level of S* or S. I. They form a group because 
they share similar characters in the shape of the found cluster and epidermal cells, 
vesture on stipe bases, and spore characters such as the perine wing and its sculp¬ 
ture. However, the simple forking pattern of the vein and brownish color of the 
spore differentiate the A. thelypterioides from A. pycnocarpon, which, in addition to 
differences in the above two characters, exhibits very peculiar simple cutting of 
lamina and pinnae. The pinnae margin of A. pycnocarpon are almost entire, with 
the veins ending in notches, not teeth. 

The closeness of these taxa is not evident from the phylogenetic diagram. 
However, it is clearly shown that A. pycnocarpon is more specialized than A. thely¬ 
pterioides, as is evidenced by the higher value of specialization index (S.I. = 8). In 
cluster analysis rather medium similarity values of .390 (T.M.) and .560 (F.M.) are 
shown between them. A. pycnocarpon has an internal cohesion value of .775 (T. M.) 
and .883 (F. M.) while A. thelypterioides has .882 (T. M.) and .936 (F. M.). They are 
both considered to be highly homogenous. 

The other group is composed of the lady fern complex. Both studies clearly 
separated the alpine lady fern from all other Ladies. In the phylogenetic scheme 
the alpine lady fern attains a high S.I. value of 8. The main differences are found 
in its skeleton like appearance, reduced condition of indusium, and structure of the 
spore. 

The point of branching is not exactly the same in both phenetic and phylogene¬ 
tic schemes. In the dendrograms the separation of alpine lady fern at a low level 
from the whole group of Ladies suggested that the former has common ancestry 
with the latter in the long past. The phylogenetic scheme, on the other hand, 
suggests that the alpine lady fern is more distantly related to the eastern American 
Ladies and shares a common ancestry with the western American Ladies. Apart 
from some obvious characters, this can be justified by the habit of the rhizome, the 
position of the young shoot in the frond cluster, and the average similarity between 
groups, being (in T. M.) .314 (to western) and .237 (to eastern: NE .223, SE .250) 
and (in F.M.) .474 (to western: NW .456, SW .491) and .375 (to eastern: NE .356, 
SE .393). Internal cohesion values are pretty high within alpine lady fern (.724 in 
T.M., and .843 in F.M.). 

This example illustrates the fact that phenetic dendrograms cannot distinguish 
between recent, rapid divergence and long-continued slow divergence. They show 
the amount of difference, but not the time of its origin. Such information will be 
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uncovered if we study and compare a series of specimens of ancient fossils and 
recent living plants. As mentioned before, such a study is not feasible because 
there is no fossil record in Athyrium. 

There remains the A. filix-femitia complex. Both schemes distinguish four 
groups. In T. M. of phenetic studies, the common ancestor gave rise to three inde¬ 
pendent groups, the northeastern, the southeastern and the western. The western 
plants divide further into two groups at a high value of S* (.60), but meaningless 
groups were formed at S* with further higher values. Hence, they are not two 
natural groups and the present author considers them to be members of the same 
subspecies — A. filix-femina subsp. cyclosorum. This is also manifested in the phy¬ 
logenetic scheme since only a portion of the spore characters (color and sculpture) 
are used for the justification of their separation. 

Within the scheme of F.M. four groups are formed (at S* = .70) directly from 
A. filix-femina complex. This means that they are very closely related, and if the 
results were taken literally, we would have to recognized all four groups at the 
same taxonomic level of hierarchy below species. As in the case of T. M. higher 
values of S* upset the schemes. Judging from the unstable nature of western 
materials in both T. M. and F. M., and rather trivial characters used for their sepa¬ 
ration in phylogenetic studies (together with reasons as outlined in the discussion 
of computer analysis), the present author considers all of the western ladies as a 
single subspecies, substantially different from both eastern subspecies. The main 
differences are in the habit of the rhizome, the presence of paraphyses in the sori, 
and the presence of stalked glands on rachis and costae. Internal cohesion shows 
that these three groups have the following high values: northeastern (T.M.: .724, 
F. M.: .838), southeastern (T.M.: .646, F.M.: .783), and western (T.M.: .628 

(SW: .785, NW: .755), F.M.: .869 (SW: .879, NW: .860)). 

The separation of eastern ladies is quite clear cut. As is evident from the 
phylogenetic diagram, they are more closely related to each other than to the 
western materials, although this is not depicted in the phenetic dendrograms. 

The clustering techniques we used here provide an aid, just like morphological, 
cytological and biochemical data, to the taxonomist, but they are not a substitute 
for the specialist. As has been shown in the above discussion, it does not make 
decisions about the taxonomic level of any group. No definitions of hierarchial 
categories ar given, because there are no known rules by which one may consis¬ 
tently apply such definitions to clusters of biological objects. However, in the 
present case, the output information is useful in making decisions about the nature, 
and even the hierarchy, of the taxa included in the present study. 

We know very well that the computer is nothing more than an instrument, just 
like an electron microscope. Profitable use of this hardware depends solely on the 
training and skill of a software programmer and taxonomist, either both existing 
as separate individuals or both in one. The computer can order and detect multiple 
correlations strictly according to the rules and regulations set up by the program¬ 
mer himself. But it can do the job much faster and more accurately that the 
human mind, and in a precisely repeatable fashion. On the other hand, it is expen¬ 
sive and time consuming and requires subjective taxonomic judgment and inter¬ 
pretation on data input and output. When taximetric methods are coupled with the 
philosophy of evolution it will come up with better classification schemes. It is 
necessary that we should be able to discover just what thought processes are in¬ 
volved in the problem of classification and to devise means of approaching the 
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conceptualization until the sets of rules and techniques become logically consistent 
by intensely examining the validity of each (Rogers and Appan, 1969). 

To accomplish this, much revision in the training of taxonomists will have to 
be made. They have to be predisposed to both the machinery of the computer and 
its programming languages (both the hardwares and the softwares) and the princi¬ 
ples of systematics, including the vigorous logic applied in exact sciences such as 
biostatistics and mathematics. By that time, one would appreciate better the nature 
of a computer aided biosystematic work and would try to collaborate constructively 
for better systematics instead of antagonizing each other (Burtt, Hedge and Stevens, 
1970). 

By using both of these computer aided taximetric and phylogentic methods to 
synthesize and analyse comparative information about Athyrium, we not only confirm 
conclusions already reached by conventional taxonomic methods, but also greatly 
improve the reliability and expressibility of the correlation of taxa being compared. 
We believed that with studies of this kind with computer methods coupled with the 
philosophy of evolution, we have come up with classifications that are clearly more 
reliable and convincible than those based solely on one sided consideration. 

LITERATURE CITED 

ANDERSON. E.. 1941. The technique and use of mass collections in plant taxonomy. Ann. 
Missouri Bot. Gard., 28: 287-292. 

_. 1949. Introgressive Hybridization. Wiley & Sons, New York. 

BLACK1TH, R. E„ & R. A. RBYMBNT, 1971. Multivariate Morphometries. Academic Press, Lon¬ 
don. 

BIR, S.S., 1965. Cytomorphological studies in the families Aspleniaceae, Blechnaceae, and 
Vittariaceae. Ph. D. Thesis, Punjab Univ, India. 

BURTT, B.L., I.C. HBDGB. & P.F. Stevens, 1970. A taxonomic critique of recent numerical 
studies in Ericales and Salvia. Notes Roy. Bot. Gard. Edinb., 30: 141-158. 
butters, F. K., 1917. Taxonomic and geographic studies in North American Ferns. I. The 
genus Athyrium and the North American ferns allied to Athyrium filix-femina. Rhodora, 
19: 169-216. 

CAM1N, J. H, & R. R. SOKAL, 1965. A method for deducing branching sequences in phylogeny. 
Evolution, 19: 311-326. 

CAVAUI-SFORZA, L.L., & A.W.F. EDWARDS. 1967. Phylogenetic analysis: models and estima¬ 
tion procedures. Evolution, 21: 550-570. 

CLUTE, W.N., 1938. Our Ferns: Their Haunts. Habits and Folklores. 2nd Ed. Stokes Co., New 
York. 

CRONQUIST, a., 1968. The Evolution and Classification of Flowering Plants. Mifflin Co, New 
York. 

DANSER, B.H, 1950. A theory of systematics. Bibl. Biotheor, 4: 115-180. 

FARRIS. J.S.. 1970. Methods for computing Wagner Trees. Syst. Zool.. 19 : 83-92. 

FERNALD. M. L„ 1946. Some trivial American forms of lady fern. Rhodora, 48 : 389-91. 
FITCH, W. M, & E. MARGOLIASH, 1967. Construction of phylogenetic trees. Sciences, 155: 
279-284. 

_, 1970. The usefulness of amino acids and nucleotide sequences in evolu¬ 
tionary studies. In Th. Dobzansky. M. K. Hecht. & W. C. Steere (eds.) Evoi. Biol., 4: 67-109. 
Meredith Pub. Co.. New York. 

HARDIN. J.W., 1937. A revision of the American Hippocastanceae. Brittonia, 9: 145-171. 
HAUKE, R. L, 1963. A taxonomic monograph of the genus Equisetum subgenus Hippod.aete. J. 
Cramer, Weinheim. 

HOLTTUM, R. E, 1940. The fern genus Diplatium in the Malay Peninsula. Gardens Bull, 11: 



Oct., 1972 Liew—Studies on North American Lady Ferns and their Allies 215 


HORNE, S.L.. 1967. Comparison of primate catalase tryptic peptides and implications for the 
study of molecular evolution. Evolution, 21: 771-786. 

ILTIS. H. H., 1959. Studies in the Capparidaceae. VI. Cleome sect. Physostenum: Taxonomy. 

geography, and evolution. Brittonia, 11: 123-162. 

JARDINE, N„ & R. SIBSON. 1971. Mathematical Taxonomy. Wiley & Sons, London. 

KAYE, R„ 1965. Variation in Alhyrium in the British Isles. Brit. Fern Gaz., 9: 197-204. 
Keener, C.S., 1967. A biosystematic study of Clematis Subsection Integri/oliae. J. Elisha Mit¬ 
chell Soc., 83: 1-42. 

KLUGE, A.G., & J.S. FARRIS, 1969. Quantitative Phyletics and the evolution of anurans. Syst. 
Zool., 18: 1-32. 

LBGENDRB, P.. & P. Vaillancourt, 1969. A mathetlcal model for the entities species & genus. 
Taxon, 18 : 245-252. 

LOCKHART. W.R., & J. Liston, 1970. Methods for Numerical Taxonomy. Vols. 1 to 4. Ameri¬ 
can Society for Microbiology, Bethesda, Maryland. 

LOWE. E.J., 1876. Our Native Ferns. Bell & Sons. London. 

MEHRA, P. N„ & S. C. VBRMA. 1954. The Cytology of some Alhyrium species from Northern 
India. Ann. Bot., n. s., 21: 455-459. 

MIKEL, J.T., 1982. A monographic study of the fern genus Anemia, subgenus Coptophyllum. 
Iowa State J. Sci.. 36 : 349-482. 

MORTON. C. V., 1963a. In H. A. Gleason's (ed.) New Britton & Brown Illustrated Flora of the 
Northeastern U.S. & Adjacent Canada. Vol. 1, p. 42-43. New York Botanical Garden, New 
York. 

_ 1963b. In H. A. Gleason & A. Cronquist’s Manual of Vascular Plants of Northeastern 

U.S. & Adjacent Canada, p. 19. Hafner Pub. Co., New York. 

PRANCE. G. T., 1967. In J. Rubin's "Optimal classification into groups an approach for solving 
the taxonomy problem.” J. Theoret. Biol., 15: 103. 

_D.J. ROGERS, & F. WHITE. 1969. A taximetric study of an Angiosperm family: 

Generic delimitation in the Chrysobalanaceae. New Phytol., 68: 1203-1234. 

ROGERS, D.J. & T.T. TANIMOTO. 1960. A computer program for classifying plants. Science, 
132: 1115-1118. 

_, H. S. FLEMING, & G F. ESTABROOK, 1967. Use of computers in studies of taxonomy 

& Evolution. In Th. Dobzansky. M.K. Hecht & W.C. Steere (eds.) Evol. Biol., 1: 169-196. 
Meredith Pub. Co.. N. Y. 

ROGERS, D.J., & S.G. APPAN, 1969. Taximetric methods for delimiting biological species. 
Taxon. 18: 609-752. 

RUBIN, ]., 1966. An approach to organizing data into homogenous groups. Syst. Zool., 15: 169. 
RUBIN, J.. 1967. Optimal classification into groups: An approach for solving the taxonomy 
problem. J. Theoret. Biol., 15: 103. 

SHAVER, J. M.. 1954. Ferns of Tenneessee. Geo. Peabody College for Teachers, Nashville, Tenn. 
SLEDGE, W. A., 1962. The Athyrioid ferns of Ceylm. Bull. British Mus. (Nat. Hist.) Bot, 2: 
277-323. 

SMALL. J. K„ 1938. Ferns of the Southeastern States. Science Press Printing Co, Lancaster, 
Pa. 

SNEATH. P. H. A, 1957. Application of computers to taxonomy. J. Gen. Microbiol, 17: 201-226. 

_ & R. R. SOKAL. 1962. Numerical taxonomy. Nature. 193 : 855-860. 

SOKAL, R.R, & P.H. A. SNEATH. 1963. Principles of Numerical Taxonomy. Freeman Book 
Co, San Francisco. Calif. 

SOLBRIG, O.T, 1970. The phylogeny of Gutkrrezia: An eclectic approach. Brittonia, 22: 217- 

STERN, K.R.. 1961. Revision of Dicentra (Fumariaceae). Brittonia. 13: 1-57. 

TRYON. R.M.. Jr, N.C. FAS3ETT, D.W. DUNLOP. & M.E. DIEMAR, 1953. The Ferns & Fern 
Allies of Wisconsin. 2nd ed. Univ. Wisconsin Press. Madison, Wisconsin. 

Wagner, W. H. Jr, 1951. Problems in the classification of ferns. In Recent Adv. in Bot, 1: 



TA1WANIA 


:>1. 17, No. 2 


841-844. Univ. Toronto Press., Toronto. 

_. 1969. The construction of a classification. In Systematic Biology, p. 67-90. Nat 

Acad. Sci., Pub. No. 1692. 

Weatkerby, C.A.. 1936. A list of varieties & forms of the ferns of eastern North America. 
Amer. Fern J„ 26: 130-136. 

Wherry, E.T., 1948. Remarks on the American Lady Ferns. Amer. Fern J., 38: 155-158. 

_1961. The Fern Guide. Doubleday & Co., Garden City, New York. 

WHIFFIN, T., & M. W. BIERNER, 1972. A quick method for computing Wagner trees. Taxon, 
21: 83-90. 

WILSON E. O., 1965. A consistency test for phytogenies based on contemporaneous species. Syst. 
Zool., 14: 214-220. 

WIRTH, M., G. ESTABROOK, & D. J. ROGERS, 1966. A graph theory model for systematic biology 
with an example for Oncidiineae (Orchidoceae). Syst. Zool., 15 : 59-69. 


APPENDIX I. List of 99 characters and their states. 







TAIWANIA 


Vol. 17, No. 2 






Oct, 1972 


Liew—Studies on North American Lady Ferns a 


Allies 


219 


DL Deltoid-lanceolate 
BL Oblong-lanceolate 
EL Elliptic-lanceolate 
EO Elliptic-oblong 
OB Oblong 
LO Linear-oblong 
EE Elliptic 

37. Segments: Shape of segment tip (SGT1P) 

LA Long-acuminate 
LU Long-acute 
AM Acuminate 
AC Acute 
OB Obtuse 
TR Truncate 
AP Apiculate 
CA Caudate 

38. Segments: Margin (SGMARG) 

EN Entire 
NO Notched 
TO Toothed 
PF Pinnatifid 
PN Pinnate 

39. Segments: Base (SGBASE) 

BI Broadly inserted 
CO Constricted 
SC Strongly constricted 
PE Pedicelled 

40. Segments: Density (SGDENSE) 

CL Close 
DI Distant 

41. Utimate segments: Shape (USGSHAP) 

LI Linear 

LL Linear-lanceolate 
LA Lanceolate 
OL Ovate-lanceolate 
DL Deltoid-lanceolate 
BL Oblong-lanceolate 
EL Elliptic-lanceolate 
EO Elliptic-oblong 
OB Oblong 
LO Linear-oblong 
EE Elliptic 

42. Ultimate segments: Shape of tip 

(USGTIP) 

LA Long-acuminate 
LU Long-acute 
AM Acuminate 
AC Acute 
OB Obtuse 

AP Apiculate 
CA Caudate 
'13. Ultimate segments 
EN Entire 


NO Notched 
TO Toothed 
PF Pinnatifid 
PN Pinnate 

44. Ultimate segments: Base (USGBASE) 

TR Truncate 
CU Cuneate 
CO Cordate 
LO Lobed 

45. Ultimate segments: Density (USGDENSE) 

CL Close 
DI Distant 

46. Cobwebby long hair, in young shoot 

petiole base (COB-HAIR) 

NO None 
FE Few 
AB Abundant 

47. Cobwebby long hair in petiole base 

(COBHA-PB) 

NO None 
FE Few 
AB Abundant 

48. Cobwebby long hair in petiole and rachis 

(COBHA-PR) 

NO None 
FE Few 
AB Abundant 

49. Short brown hairs, in groove of petiole 

(SBH-PG) 

NO None 
FE Few 
AB Abundant 

50. Short brown hairs in groove of rachis 

(SBH-RG) 

NO None 
FE Few 
AB Abundant 

51. Short brown hair in pinnae rachis 

(SBH-PR) 

NO None 
FE Few 
AB Abundant 

52. Slender long hairs, on pinnae rachis 

(SLH-PR) 

NO None 
FE Few 
AB Abundant 

53. Star shaped short hairs, on rachis 

(SSH-RA) 

NO None 
FE Few 
AB Abundant 

54. Star shaped short hairs in pinnae axils 

(SS1I-PA) 


): Margin (USGMAKG) 
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NO None 
FE Few 
AB Abundant 

55. Star shaped short hairs on pinnae rachis 

(SSH-PR) 

NO None 
FE Few 
AB Abundant 

56. Short glandular hair on undersurface of 

blade (SGHBLAD) 

PR Present 
AB Absent 

57. Density of all kinds of hairs, on petiole 

bases (HAIR-PB) 

NO None 
FE Few 
AB Abundant 

58. Density of all kinds of hairs on rachis 

(HAIR-RA) 

NO None 
FE Few 


59. Density of all kinds of hairs on pinnae 

rachis (HA1R-PR) 

NO None 
FE Few 
AB Abundant 

60. Density of all kinds of hairs in pinnae 

axil (HAIR-PA) 

NO None 
FE Few 
AB Abundant 

61. Density of all kinds of hairs on costae 

(HAIR-COS) 

NO None 
FE Few 
AB Abundant 

62. Density of all kinds of hairs on the leaf 

blade as a whole as seen by eye 
(HA-BLAD) 

NO None 
FE Few 
AB Abundant 

63. Wide ovate thin scales found on petiole 
base (WOTS-PB) 

PR Present 
AB Absent 

64. Broad and short scales found on petiole 

base (BSS-PB) 

PR Present 
AB Absent 

65. Long slender scales on petiole base, 

rachis, costae, pinnae axil 
(LSS-PRCA) 


PR Present 
AB Absent 

5. Long slender scales: Size—length (mm) 
(LSS-LENG) 

1. Long slender scales: Size—width (mm) 
(LSS-WIDE) 

i. Long slender scales: Structure of scales 
(LSS-STRU) 

FI Fibrous 
TW Thin-walled 

). Long slender scales: Color of scales 
(LSS-COLO) 

YE Yellowish concolorous 
YB Yellowish brown 
BB Brownish black 

). Long slender scales: Color of scale at 
center (LSS-CECO) 

BL Blackish 
DB Dark brown 

CO Concolorous with other part 
l. Long slender scales: 

Lengt 


Widtl 


2. Slender and sh 
pinnae axis. 

(SSS-RPCA) 

PR Present 
AB Absent 

I. Density of all kinds of scales 
base (SCALE-PB) 

NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 

I. Density of all kinds of scales 
(SCALE-RA) 

NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 

i. Density of all kinds of scales: 
rachis (SCALE-PR) 

NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 

i. Density of all kinds of scales: 
axil (SCALE-PA) 

NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 
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NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 

78. Density of all kinds of scales: on the 

(SCA-BLAD) 

NO Glabrous 
SC Scanty 
MA Many 
AB Abundant 

79. Shape of sori (and indusium): asplenioid 

sori (SORI-ASP) in % 

89. Shape of sori (and indusium): athyrioid 
(SORI-ATH) in % 

81. Shape of sori (and indusium): diplazioid 

(SORI-DIP) in % 

82. Shape of sori (and indusium): roundish 

(SORI-ROU) in % 

83. Size of asplenioid sori (and indusium): 

Length (ASP-LENG) in mm 
81. Size of asplenioid sori (and indusium): 
Width (ASP-WIDE) in mm 

85. Diameter of roundish sori (RSORIDI) in 

86. Position of sori on the vein (SORIPOSI) 

MA Marginal 
SU Submarginal 
ME Medial 

OV Along oblique vein 

87. Transitional sequence of sori 

(SORISEQU) 

SO Asp—»dip—>ath—»asp 
ST Ath—»dip—»asp 
NO None series formation 

88. Indusium presence (INDUSIUM) 

PR Present 
RE Remant 
AB Absent 

89. Free margin of indusia (INDUMARG) 

EN Entire 

ET Entire with a few irregular teeth 
DE Dentate 
SC Short cilia 
LC Long multicellular cilia 
LG Long cilia with glands 

93. Length of indusia cilia (INDUCILI) 

SH Short 


LO Long 
VL Very long 

91. Color of the glands of indusia cilia 

(INCICOLO) 

YB Yellowish brown 
RD Red 
GR Green 

92. Length of sporangial sulk (at 90 x) 

(STAKLENG) number of grid 
9a Structure of sporangial stalk 
(SPORSTAK) 

NP Non proliferate 

PR Proliferate with no gland 

PG Proliferate with glandular structure 

94. Spore color (SPOCOLOR) 

YE Yellowish 
YB Yellowish brown 
BR Brown 
BL Black 

95. Spore size-Length (p) (SPOLENG) 

96. Spore size-Width (p) (SPOWIDE) 

97. Exospore wing (SPOWING) 

VL Very long 
PR Prominent 
BO Barely observable 
AB Absent or hardly observable 

98. Perine sculpture (SPOSCULP) 

SM Smooth 
SP Slightly papillate 
HW Heavily warty 
RE Reticulate 
RI Ridged 
FO Folded 

99. Range (RANGE) 

EA Eurasia 

NE Northeastern North America 
SE Southeastern North America 
NW Northwestern North America 
SW Southwestern North America 
EE Northeastern and southeastern 
North America 
UE Eurasia, northeastern and 
southeastern North America 
AM Northeastern, northwestern and 
southwestern North America 
WW Northwestern and southwestern 
North America 





